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ABSTRACT 

Fluoride  concentrations  in  complex  natural  waters,  such  as  estu¬ 
aries,  rivers,  or  lakes,  may  be  measured  using  an  ion-selective 
electrode  and  a  known-addition  technique  in  conjunction  with  a  com¬ 
pleting  buffer.  No  additional  information,  such  as  interfering  ion 
concentrations  or  ionic  strength,  is  needed.  Reliable  measurements 
may  be  made  in  river  waters  containing  as  much  as  25  mg/1  Al3*, 
125  mg/1  Fe3t,  20  g/l  Ca24,  or  2.5  g/l  Mg2*.  Results  are  accurate 
to  about  0.02  mg  F/l.  In  rivers,  the  relative  standard  deviation  is 
1%  to  2%  when  the  fluoride  concentration  is  1  mg/1  and  up  to  5%  when 
its  concentration  is  0.1  mg/1.  In  estuaries,  the  relative  standard 
deviation  ranges  between  0.3%  and  1,0%  under  typical  conditions. 
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DETERMINATION  OF  FLUORIDE  POLLUTANTS  IN  NATURAL  WATERS 
USING  A  KNOWN  ADDITION  TECHNIQUE 


INTRODUCTION 

Fluorine  is  an  mportant  toxicant  that  must  be  monitored  in  the  environment.  A 
fluoride  concentration  of  more  than  1.  2  to  1.  5  mg/1  in  drinking  water  causes  damage, 
while  too  little  promotes  dental  caries.  Fluorides  are  classified  as  highly  toxic  (Sax, 
1983);  they  harm  both  plants  and  animals  when  introduced  into  the  atmosphere  as  gases, 
mainly  HF  and  SiF4,  and  particulates  or  when  injected  into  surface  waters  (McCune, 

1969;  Hunter,  1962;  Shupe,  1969).  High  local  concentrations  can  come  from  volcanic 
emanations,  from  industrial  plants  processing  phosphate/rock/fertilizer,  steel,  alumi¬ 
num,  chemicals,  coai/fuei  oil,  brick/tile/ceramics,  and  uranium  (Carpenter,  1969; 
Hunter,  1962) ,  and  from  feed  supplements,  insecticides,  and  rodenticides  rich  in  fluo¬ 
rides  (Shupe,  1069).  Most  local  emissions  of  airborne  fluoride  quickly  cause  increased 
levels  in  nearby  vegetation  or  surface  waters,  which  then  affect  animal  life.  In  one 
instance,  surface  waters  10  km  from  an  industrial  airborne  source  contained  fluoride 
concentrations  of  10.9  mg/1,  and  within  a  5-km  distance,  trees,  vegetables,  fruit, 
bee  colonies,  and  cattle  were  adversely  affected,  and  children’s  blood  showed  decreased 
hemoglobin  and  increased  erythrocyte  levels  (Marier,  1968). 

To  assess  the  extent  and  effects  of  fluorine  pollution,  measurements  are  needed 
in  a  wide  range  of  media,  such  as  air,  water,  plant  and  animal  tissues,  bones,  minerals, 
and  body  fluids  such  as  urine  and  blood.  Fluoride  analysis  in  complex  solutions  has 
been  rather  difficult  until  recently;  the  usual  procedures  required  considerable  manipu¬ 
lation  to  separate  fluoride  from  interfering  substances.  Simpler  techniques  became 
possible  (Crosby,  1968)  with  the  recent  invention  of  fluoride -selective  electrodes,  for 
such  electrodes  respond  directly  to  free  fluoride  ion  activity,  even  in  the  presence  of 
large  excesses  of  most  other  ions.  They  have  been  used  for  the  determination  of  fluo-ride 
in  a  variety  of  materials,  such  as  phosphate  rocks  (Edmond,  196S),  bones  (Singer  and 
Armstrong,  1968),  air  and  stack  gas  (Elfers  and  Decker,  1968),  seawater  (Warner, 

1969a;  Brewer,  Spencer  and  Wilkniss,  1970),  and  potable  waters  (Frant  and  Ross, 

1968;  Light,  Mannion,  and  Fletcher  1969). 

Most  measurements  have  used  direct  potentiometry.  It  is  rapid  and  convenient, 
but  the  solution  chemistry  must  be  known  or  controlled  if  the  total  concentration  is  to  be 
determined,  since  the  electrode  responds  only  to  free  ion  activity.  Any  fluoride  bound 
in  complexes  will  not  be  detected,  and  the  ionic  strength  of  the  standard  must  match 
that  of  the  unknown  so  that  the  activity  coefficients  are  the  same.  Many  ingenious 
techniques  have  been  developed  in  which  the  chemical  compositions  of  samples  and 
standards  are  closely  matched  or  in  which  high-ionic-strength  complexing  buffers  are 
used  to  free  complexed  fluoride  and  to  fix  the  ionic  strength  of  the  sample  at  some 
known  high  value.  Harwood  (1969)  showed  that  the  buffer  suggested  by  Frant  and  Ross 
(1968)  was  effective  in  low-ionic-strength  river  water  if  L Ai+3J  was  low,  but  that  signi¬ 
ficant  interference  occurred  if  aluminum  was  present  at  concentrations  of  0. 2/mg/l  when 
the  fluoride  concentration  was  1  mg/1.  Using  an  improved  buffer,  with  cyclohexanedia- 
mine  tetraacetic  acid  as  the  trivalent  ion  complexing  agent,  he  obtained  95%  recovery  in 
water  containing  3  mg  Al/1. 
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The  i  .‘nerality  of  the  complexing  buffer  approach  is  limited,  however.  Actual 
,  A  l’3  m*.  l  he  known  (or  be  known  to  be  small)  to  assess  the  probable  errors  of  potentio- 
metric  n  ;asurenients,  and  these  errors  will  become  more  important  in  untreated 
rivers,  There  F  is  lower.  These  techniques  become  increasingly  difficult  to  apply 
when  co  mlexing  ions  are  present  in  high  and  unknown  amounts,  or  when  the  total  ionic 
strengt  pi  the  sample  is  appreciable.  This  is  the  situation  in  some  polluted  rivers  and 
lakes,  ,nd  particularly  in  estuaries  and  saline  lakes. 

T)  e  main  ions  in  natural  waters  that  complex  fluoride  ions  are  Ca42,  Mg+2,  Al43, 
and  Ft  ’3.  The  divalent  ions  form  weak  complexes  and  the  trivalent  ions  form  very 
strom.  ones.  In  estuaries,  the  water  is  essentially  seawater  as  diluted  by  river  water 
(Carpenter,  1957).  Magnesium  is  abundant,  iron,  aluminum,  and  calcium  are  low,  ,F_] 
is  ab  >ut  0.  2  to  1.4  mg/1,  and  ionic  strengths  range  up  to  0.7.  In  rivers,  calcium 
predf  minutes  over  magnesium,  and  variable  amounts  of  iron  and  aluminum  are  found. 
Typi  fal  values  of  iron  and  aluminum  concentrations  in  rivers  are  0. 1  mg/1,  and  the 
high  ist  concentrations  encountered  are  unlikely  to  exceed  25  mg/l  Fe  and  10  mg/1  Al. 
The  concentration  of  dissolved  F  is  about  0. 1  to  0.  2  mg/l,  although  certain  exceptional 
values  as  high  as  23  mg/l  have  been  reported  (Livingstone,  1963).  Very  wide  variations 
are  found  in  saline  lakes. 

In  such  solutions,  considerable  simplification  can  result  by  using  a  known-addition 
or  spike  technique  (Garrels,  1967;  Durst,  1969;  Orion,  1969;  Orion,  1970),  and  in  some 
ca  ses  it  may  be  the  only  feasible  approach.  In  the  known-addition  method,  the  unknown 
simple  itself  becomes  the  chemical  matrix  for  the  standardizing  sample.  The  ionic 
e ^engtti  can  be  very  high,  and  its  value  need  not  be  known.  Very  high  concentrations 
(  t  i  ns  that  complex  fluoride  can  be  tolerated,  and  in  most  cases  their  identities  need 
oe  known.  Accurate  determinations  are  possible  even  when  the  majority  of  the 
uur.de  is  bound  in  complexes.  Although  Bock  and  Strecker  (1968)  and  Baumann  (1968/ 
/if  monstrated  the  usefulness  of  spike  methods  in  waters  of  different  ionic  strengths  and 
I  c  <  taming  various  interfering  ions,  these  methods  have  been  little  used  to  date  in  com- 
pu  x  natural  fluids.  The  purpose  of  this  study  was  to  devise  a  method  of  analysis  based 

/  n  known  additions  that  would  be  generally  useful  down  to  0. 1  mg  F/l  in  various  natural 
waters,  with  no  prior  knowledge  of  chemical  composition.  The  principles  involved  are 
i  general  and  apply  to  other  ion-selective  electrodes  as  well,  and  while  this  application 
)  is  to  natural  waters,  analogous  techniques  should  prove  useful  for  a  wide  range  of 
aqueous  solutions  prepared  from  a  number  of  sources. 


THEORETICAL  CONSIDERATIONS 

If  c0  is  the  total  concentration  of  the  ion  of  interest  in  the  unknown  and  CA  is  the 
amount  its  concentration  is  increased  by  the  spike,  then  according  to  the  formulation  of 
Orion  (1969,  1970), 

C0  CA  (eAE/S  t)  (1) 

where  -\E  is  the  potential  change  observed  on  spiking  and  s,  the  response  slope  of  the 
electrode,  is  the  theoretically  predicted  RT/F  for  the  fluoride  electrode  (Warner, 

1969b).  This  equation  assumes  that  (a)  the  change  in  total  ionic  strength upon  addition 
of  the  spike  can  be  neglected,  (b)  the  fraction  of  the  ion  being  measured  that  is  free  and 
uncomplexed  remains  essentially  constant,  and  (c)  electrode  interferences  are  not  pres¬ 
ent  in  amounts  that  will  affect  electrixie  response.  If  the  concentration  of  the  spike 


/ 


NRL  REPORT  7216 


3 


solution  is  Ci  and  the  volume  of  test  solution  and  added  spike  are  v0  and  V),  respec¬ 
tively,  then 


Ci  vi  ~  co  V' 
v0  +  V 


or,  to  a  good  approximation  when  c,  >  100  c0, 


CA  = 


C.  vt 

V0  +  V 


(3) 


In  natural  waters,  assumptions  (a)  and  (c)  are  fulfilled  by  using  a  high-ionic-strength 
buffer,  keeping  Vi  <  0.01  v0  and  Ca«  C0.  In  most  solutions,  assumption  (b)  is  also  justi¬ 
fied.  As  is  shown  below,  the  degree  to  which  it  is  justified  depends  on  the  amounts  and 
kinds  of  complexing  ions  present. 

Let  p  be  the  fraction  of  the  total  fluoride  present  in  free,  uncomplexed  form  when 
in  the  presence  of  complexing  ion  Mn+,  whose  first  concentration  equilibrium  constant 

K,.  is 


K  = 

iMn+]  m 

Since 

t,  — jo _ 

If']  +  ] 


then  combining  Eqs.  (4)  and  (5)  gives 


P  = 


1 


1  +  K„  ( M' 


(4) 


(5) 


(6) 


When  F  ions  are  added  in  the  spike,  some  must  combine  with  M"4  ions,  causing  ; 
to  increase.  If  4  increases  by  1%,  c0  calculated  by  Eq.  (I)  will  be  2%  low  under  the 
stated  conditions  (Appendix  A).  Taking  this  value  as  an  upper  limit  for  error  tolerable 
from  a  change  in  <t>,  solution  conditions  may  be  defined  under  which  a  spike  technique 
can  be  used.  H  KclMn4]  <  0.01,  then  any  decreases  in  lMn‘]  upon  spiking  cannot  change  ■; 
more  than  1%.  If  Kc  (Mn‘]  >  o.oi ,  two  cases  can  arise.  If  [Mn*]  then  the  addition 

of  F"  will  cause  small  changes  intMn*l  and  again  bp  <  1%.  However,  if  iMn4!«  if  ]  «  to"  m 
(for  natural  waters)  &nd  if  ->o.oi,  then  A^can  exceed  1%,  but  this  requires  that 

Kc  ■  ioJ.  For  the  fluoride  case,  then,  large  background  concentrations  of  ions  such 
as  Mg2*(Kt  -  18  for  „  0.5)  and  Ca2’(K(.  3  for  «  0.5)  will  not  cause  appreciable  error. 

The  only  ions  that  can  interfere  are  the  relatively  few  with  Kc  io3,  such  as  Al3*  and 
Fe3'.  These  ions  can  be  reduced  to  sufficiently  low  concentrations  by  using  the  com¬ 
plexing  reage  '  suggested  by  Harwood  (1969).  In  questionable  cases,  results  may  be 
checked  using  a  double  spike  (Orion,  1970). 
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EXPERIMENTAL  PROCEDURE 
Apparatus 

Potentials  between  an  Orion  94-09  LaF3  electrode  and  a  saturated  calomel  electrode, 
thermostated  at  25.0  ±  0. 1°C,  were  measured  with  an  Orion  Model  801  pH  meter. 
Filtered  voltage  was  displayed  with  suppressed  zero  on  a  recorder,  5  mV  full  scale. 
Spikes  ranging  from  0. 05  to  0.  25  ml  were  added  using  a  Manostat  Digipet  accurate  to 
±  0. 001  ml. 


Reagents 

Fluoride  standards  were  prepared  by  weight  and  diluted  as  necessary.  The  complex- 
ing  buffer  (TISAB  IV)  was  a  slight  modification  of  that  proposed  by  Harwood  (1969).  It 
was  identical  to  the  TISAB  of  Frant  and  Ross  (1968)  except  that  the  sodium  citrate  was 
replaced  by  5  g  of  1,  2-cyclohexanediamine  tetraacetic  acid  (CDTA)  per  liter  of  solution. 
Spike  solutions  were  10  2M  NaF  for  estuarine  vate”  and  river  water  containing  50  n  M  F 
(1  mg/l)  and  10~3M  NaF  for  river  water  containing  5mM  F  . 

Mg(NC>3) 2,  CaCl2,  FeCl3,  and  A1 2(804)3  •  K2SO4  •  24H2  0  were  added  to  aqueous 
NaF  standards  to  make  synthetic  river  waters.  All  reagents  were  fluoride-free  except 
CaCl2,  whose  fluoride  content  was  determined  by  both  the  previously  described  method 
(Warner,  1969b)  and  by  the  spike  technique  discussed  here.  Simulated  estuarine  waters 
were  similarly  prepared  in  diluted  seawaters  whose  fluoride  concentrations  were  first 
determined  using  the  potentiometric  method  previously  described  (Warner,  1969a). 
Because  complexation  reactions  can  be  slow,  waters  were  prepared  at  least  24  hours 
before  use,  and  control  samples  were  analyzed  after  storage  up  to  2  months.  Synthetic 
river  waters  containing  high  iron  concentrations  were  made  sufficiently  acidic  to  hold 
all  iron  in  solution  prior  to  TISAB  IV  addition;  otherwise  colloidal  precipitates  of  iron 
hydroxide  irreversibly  removed  fluoride  from  solution.  In  estuarine  waters,  the  pH 
was  not  modified  and  iron  was  permitted  to  precipitate  partially  when  added  to  simulate 
iron- rich  water  pouring  into  an  estuary. 


Method 

Samples  were  mixed  with  TISAB  IV  in  the  ratio  of  5  ml  of  sample  to  1  ml  of  buffer, 
making  available  2.9  millimoles  of  CDTA  per  liter  of  test  sample.  Complexation  of 
interfering  ions  was  rapid  for  all  ions  except  A 1  * for  A 1*'  the  fraction  complexed  at 
25UC  was  still  changing  after  2  hours,  but  was  constant  after  20  hours.  Accordingly, 
all  samples  were  mixed  with  buffer  at  least  20  hours  before  analysis.  Usually,  25  ml 
of  the  mixtures  was  taken  for  analysis.  Potentials  before  and  after  spiking  (C  •c0)  were 
measured  after  15  minutes.  If  a  steady-state  potential  was  not  achieved,  the  remaining 
increment  was  estimated  from  the  voltage-time  curve;  estimates  did  not  exceed  0.3  mV. 
Using  Eq.  (1),  the  total  concentration  of  fluoride  C0  was  computed  from  the  known  and 
c  ,  using  the  values  of  c0  c  in  Table  1. 

For  5,  M  F  (0. 1  mg/l)  in  synthetic  river  water,  slow  electrode  response  made 
measurements  tedious  at  E  l«  mV.  Smaller  spikes  were  used  (  F.  4  «  mV),  ant.  a  tech¬ 
nique  sacrificing  accuracy  for  speed  was  adopted.  Potential  measurements  were  arbi¬ 
trarily  made  15  minutes  after  each  change  in  solution  composition,  and  the  remaining 
increment  to  steady-state  potential  was  estimated.  Estimates  were  limited  to  0.3  mV  to 
avoid  the  introduction  of  personal  bus,  although  true  increments  were  always  larger. 

This  procedure  consistently  underestimated  'E  and  all  resulting  values  for  C0  were 
positively  biased,  but  within  acceptable  limits.  To  a  first  approximation,  the  error 


Table  1 

Values  of  c0  c  as  a  Function  of  ,E  at  25 °C  for  Electrodes 
with  Theoretical  Monovalent  Slope  of  S  =  RT  F* r 
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■  Corresponding:  values  of  Q,  C  for  electrodes  with  theoretical  divalent  slcpe  (S  -  RT  2F)  can  be  touna  by  doubling  tne  odsc.  ea 
using  the  above  table;  similar  proportioning  allows  its  use  with  electrodes  having  experimentally  determined  slopes,  or  at  other 
temperatures. 

The  value  of  E  is  the  number  in  this  column  plus  the  number  at  the  top  of  any  particular  column  to  the  right. 
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In  c0  will  be  about  equal  to  the  negative  percent  error  in  E;  e.  g. ,  if  E  is  lO'T  low,  c 
will  be  about  10%  high. 


RESULTS  AND  DISCUSSION 

Results  of  fluoride  analyses  obtained  in  seawater  and  simulated  estuarine  waters 
ranging  in  salinity  from  7  to  35  %>o  are  summarized  in  Table  2.  Concentrations  <t  in¬ 
terfering  ions  were  selected  well  above  those  expected  from  mixing  river  waters  with 
estuary  waters.  Because  Mg2’ complexes  much  more  F  than  Ca2',  the  former  was  the 
main  divalent  cation  interference  studied.  No  difference  was  observed  for  measure¬ 
ments  made  in  synthetic  and  real  seawaters,  and  only  the  results  for  real  seawaters 
are  reported. 

Results  in  simulated  river  waters  are  shown  in  Tables  3  and  4.  The  concentrations 
of  interfering  ions  studied  were  well  above  those  normally  encountered  in  natural  waters, 
even  when  severely  polluted.  Since  CDTA  was  only  available  to  complex  2.9  millinmU  s 
per  liter  of  trivalent  ions,  concentrations  above  that  level  resulted  tn  no  measurable  F 
It  is  probable  that  higher  concentrations  of  CDTA  would  permit  fluoride  determinations 
in  the  presence  of  even  larger  amounts  ol  interfering  ions,  but  this  has  not  been  experi¬ 
mentally  verified.  Results  in  Tabic  3  are  representative  of  high-fluoride  waters  or 
those  in  which  fluoride  has  been  increased  by  trea!  nsnt.  Results  in  Table  4  are  typical 
for  untreated  rivers.  The  actual  error  to  be  expected  is  nearly  the  sar.e  m  each  water. 
A  2%  error  at  the  50  mM  level  is  1  uM/l  as  is  a  20%  error  at  5  ..M/1. 

In  the  high-fluoride  waters,  a  significant  error  (of  about  8%)  is  noted  only  for  waters 
containing  very  high  (0.  5M)  magnesium  levels.  In  general,  at  50  ..M  F  in  virtually  any 
natural  water,  results  should  be  accurate  to  2%,  end  the  relative  standard  deviation  of  at 
individual  measurement  will  re  1  to  2%.  At  the  5.. M  F  level,  using  the  rapid  method 
described,  results  will  be  5%  to  30%  high,  and  the  relative  standard  deviation  will  range 
between  1%  and  5%. 


Table  3 

Fluoride  Recovery  in  Simulated  River  Waters  Containing 
50  ..M/1  (1  mg/1) 


Ion  and  Concentration  Added 

Mean  Recovery  (and  Uncertainty*)  (%) 

None 

100.7  (0.9) 

Ai  (111)  0. 4  mM  (about  10  mg/1) 

98.7  (1.6) 

1.0  mM 

101.0  (2.2) 

Fe  (in)  0. 4  mM  (about  25  mg/1) 

100.6  (0.4) 

1.0  mM 

102.2  (1.2) 

2.0  mM 

99. 1  (1.9) 

Ca  (II)  100  mM  (about  4000  mg/1) 

99.9(0.9) 

500  mM 

100.4  (0.4) 

Mg  (II)  100  mM  (about  2500  mg/l) 

103.7  (0.9) 

500  mM 

108.4  (1.3) 

Synthetic  river  lr 

101. 2  (1.7) 

Synthetic  river  2* 

101.2(1.7)  j 

•Uncertainty  la  the  relative  standard  deviation  of  an  Individual  measurement;  four  to 
eight  sample*  were  measured  at  each  concert  rail  on. 

’  Synthetic  river  1  contained  0.  I  mM  AI,  0.  1  mM  Fe,  and  25  mH  Mg. 

5  Synthetic  river  2  contained  0.  4  mM  AI.  0.  4  mM  Fe,  and  50  mM  Mg. 


8 


T.  B.  WARNEH 


Table  4 

ride  Recovery  in  Simulated  River  Water  Containing  5  M/1  (0. 1  mg/1) 


ering  Ion 
•ntracion 

Mean 

Recovery  (%) 

Rel.  Std.  Deviation 
(individual)  (%) 

Number  of 
Determinations 

Jone 

128 

4.9 

8 

Al(lII) 

125 

3.4 

4 

lM  Fe(III) 

106 

1.3 

4 

M  Ca(II) 

111 

1.0 

4 

M  Mg(II) 

114 

2.8 

4 

in-  esults  indicate  that  large  variations  in  total  ionic  strength  do  not  interfere  with 
the  accu  cy  of  such  determinations.  Very  large  concentrations  of  ions  that  form  com¬ 
plexes  v  h  low  •  ♦ability  constants  can  be  tolerated  if  a  measurable  amount  of  free  flu¬ 
oride  rt  ains  in  solutic  This  technique  offers  a  convenient  way  to  determine  the 
amount  fluoride  impurity  present  in  other  reagents,  because  commonly  the  impurity 
levels  a  highest  in  those  salts  whose  cations  lorm  flu  ride  complexes.  Similar  tech¬ 
niques  s:  juid  allow  an  accurate  determination  of  flu<v  ide  levels  in  many  other  media 
such  as  lant  and  animal  tissues  or  biological  fluids. 
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Appendix  A 


ERROR  IN  CALCULATED  VALUE  OF  FLUORIDE  CONCENTRATION  IF  THE 
FRACTION  UNCOMPLETED  AFTER  SPIKE  ADDITION  CHANGES  BY  1% 


In  the  original  solution,  the  totai  measured  potential  E  j  is 

Ej  =  E0  +  S  ln(;.C0),  (Al) 

where  c0  is  the  total  concentration  of  fluoride  present,  *  is  the  fraction  of  that  total  fluo¬ 
ride  present  in  free  uncomplexed  form,  and  E0  is  a  constant  representing  the  portion  of 
the  total  potential  due  to  reference  electrodes  and  internal  solutions  and  also  includes  the 
(unchanging)  fluoride  ion  activity  coefficient.  The  potential  E2  in  the  solution  after  spik¬ 
ing,  if  <f>  does  not  change  and  if  C&  =  c0,  is 

E2  =  E0  +  S  ln(*q,  +  4>c±)  (A2) 

and 

AE  =  Ej  -  Ej  =  S  in  2.  (A3) 

Hence,  from  Eq.  (1) 

“■=•5—=  ».000.  (A4) 

CA 


which  is  correct.  But  if  *>2,  the  fraction  actually  free  after  spiking,  is  1.01*,  then 
E2  =  E0  +  S  In  2.O2*C0,  AE  =  s  In  2-02,  and  C0  would  be  calculated  from  Eq.  (1)  as 

C0=0.98Ca:  (A5) 

that  is,  the  computed  result  would  be  2%  low. 
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river  waters  containing  as  much  as  25  mg/1  Al3+,  125  mg/1  Fe3+,  20  g/1  Ca2+  or  12  5 
g/a  Mg!+.” 


